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We present results for the exotic +_ and 1~ + bbg hybrids at f3 = 6.0. A leading order NRQCD action is used 
to generate the heavy quark propagators. The optimisation of operators is discussed and we compare our results 
to the predictions of recent potential calculations. 



1. INTRODUCTION 

In hybrid mesons the gluon string between the 
qq pair is excited and contributes to the quantum 
numbers of the particle. These mesons are pre- 
dictions of QCD which are not yet confirmed in 
experiment. More difficult to discriminate from 
quark model states than glueballs, it is partic- 
ularly important to provide accurate predictions 
from lattice studies. 

Most experimental candidates are for hybrids 
containing light quarks; however, the best system 
for a hybrid search may be cc or bb, where there is 
a large gap between the lowest state and the DD 
and BB thresholds respectively. Theoretically, 
heavy hybrids are easier to study; on the lattice, 
NRQCD can be used to simulate heavy quarks 
and the results can be compared with those from 
hybrid potentials. 

Here, we present results for bbg hybrids. The 
lattice hybrid operators are constructed from sta- 
ples of gauge links, U— n, in combinations corre- 
sponding to the and T\~ + representations of 
the cubic symmetry goup, as described in Lacock 
et al. Q. In early potential studies 0] it was 
found that this gluon excitation, corresponding 
to one unit of angular momentum about the qq 
axis, gives rise to the lowest hybrid meson states. 
Combined with qTq in a spin singlet (A^ + ) and 
triplet (T~) the hybrid meson states 

Tf~: l~~,0~ + ,l~ + ,2~ + 
Tf + : l ++ ,0 + -,l + -,2+- 
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are accessible, where the exotic states are shown 
in bold. In the static limit these states are all de- 
generate; introduction of quark motion, using the 
lowest order NRQCD action, will create a split- 
ting between the sets of states arising from the 
T± and T-f + gauge paths. Spin-spin and spin- 
orbit interactions are needed in order to obtain 
splittings between the individual states. These 
spin splittings, eg between the O ^ 1 h and 
2 h , will be ordered in the same way as for the 
non-hybrid 3 P states and are likely to be of order 
10 MeV. Since the statistical errors in our cal- 
culation are still at approximately 200 MeV we 
ignore the fine structure and use the lowest order 
NRQCD action including only the kinetic energy 
term. In this case there is no mixing between the 
hybrid and normal bb states. 

The correlators corresponding to the hybrid op- 
erators were computed on 500 quenched 16 3 x 48 
configurations at j3 = 6.0. A bare b quark mass, 
from bb spectroscopy, of 1.71a was used in the 
NRQCD action. The hybrid signal is dominated 
by noise after only a few timeslices, and multiple 
spatial and temporal sources are required. It is 
essential to optimise the operators in terms of the 
qq separation, R, depth of the staple, D, and the 
coefficient, c, and level, n, of the fuzzing of the 
gauge links. We use the results of a study of the 
hybrid potential Q on the same configurations to 
choose the optimal values of these parameters. 

2. THE HYBRID POTENTIAL 

The hybrid potential corresponding to the T\ 
representation, E u , was calculated in the stan- 
dard way using operators of the form U— n. In 
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Figure 1. 

order to find the optimal operator, the depth D 
was varied at each R and the maximal overlap 
of the operator with the ground state, A gs . , was 
sought. Table |f| presents the results as a function 
of D for the example of R — 6. A gs , is not very 
sensitive to D as long as D is small. A fairly uni- 
form optimal overlap was found for a wide range 
of R, where we investigated up to R = 9. 

Fuzzing the gauge links seemed to act as a 
smoothing function on the hybrid source and it 
was found sufficient to optimise the value of the 
spatial plaquette with respect to c and n in order 
to obtain their optimal values: roughly c — 2.3 
and n = 25. 

Using these optimal parameters for the gauge 
path, the hybrid potential shown in figure |l] was 
obtained; the ground state potential, A lg , is also 
shown for comparison. The most striking fea- 
ture of the potential is its flatness which suggests 
broad hybrid wavefunctions and densely packed 
states. To extract this information the potential 
was fitted to a functional form suggested by string 
models: 



V = V a . 
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with parameter c. The hybrid energy levels, E n i, 
and wavefunctions, u n i 1 where tp — u n {Y^ / R, 
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Figure 2. 

can then be extracted by solving the Schrodinger 
equation in the continuum. The lowest lying hy- 
brid levels are given in table || in GeV and the 
corresponding radial wavefunctions are shown in 
figure |[ The scale was set using the average of 
the lattice spacings from the bb P — S and 25 — 15 
splittings. Note that the first hybrid state is a 
P-state. The IP bb wavefunction is shown for 
comparison in the figure. 

The wavefunctions are broad, as expected. For 
a I ~ 1.8 fin lattice we can expect large finite 
volume errors for a ground state hybrid meson 
with r>. m . s . ~ 0.28 fin. Figure || suggests a lattice 
hybrid operator with a qq separation of 0.5 fm or 
~ 6 lattice units at /3 — 6.0 would minimise the 
overlap with the first excited state since the cor- 
responding wavefunction has a node at this point. 
However, a more reliable method to extract the 
ground state energy is to perform multiple expo- 
nential fits to a matrix of correlation functions 
obtained from different qq separations, where the 
different values of R give rise to very different 
overlaps with the hybrid excited states in the off- 
diagonal elements of the matrix. For example, 
from figure we see a suitable 2x2 matrix can 
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Figure 3. 
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be formed from R = 3 and R = 7. 

3. HYBRIDS WITH PROPAGATING 
QUARKS 

Considering the potential results we performed 
an initial calculation of the T± and hy- 
brids using a staple size of R = 3 and D = 2 at 
both the source and the sink. Figure || presents 
the results for the effective mass of the T± hy- 
brid obtained from 10, 000 sources. There is a 
large contribution from excited states for this op- 
erator. However, a plateau seems to begin from 
timeslice 6 or 7. One-exponential fits in the re- 
gion 6—10 are in agreement with two-exponential 
fits with t m i n = 3—5. Considering all these fits, 
we find the ground state simulation energy to be 
consistent with 1.0(1), as indicated in the figure. 
This corresponds to a splitting with the 3 <Si bb 
state of 0.55(10) in lattice units. We also find a 
splitting of 0.1 — 0.3 with the T l ~ + states, which 
confirms that the 1 h is the lowest exotic hybrid. 

Converting our results into physical units we 
obtain 11.5(5) GeV for the T x + ~ hybrid state. 
This is compared with the lower lying bb spectrum 
obtained from an 0{Mv i ) NRQCD action f§ and 
also from a potential study, accurate to the same 
order in v ^ in figure || Our prediction lies ler 
above the threshold for BB production but 1<t 
below that for BB**, which is predicted by flux 
tube models to be the dominant hybrid decay 
mode. A study by Manke et al.. |j| using dif- 



ferent operators on the same configurations finds 
consistent results with the hybrid lying la 
higher. Our results are also consistent with the 
predictions from the potential. 

The position of the hybrid state with respect 
to the two thresholds is not settled. An improved 
determination of the hybrid state is needed and 
we propose to do this by computing a matrix of 
correlators between different qq separations. At 
the moment the large statistical errors dominate 
the uncertainty in the prediction. However, the 
central value may be shifted significantly by fi- 
nite volume effects and to a lesser extent quench- 
ing and the systematic errors associated with 
NRQCD. 
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